function [x,iiIteration,err_var,time_cost]=Simplified_Newton_Equation_Solver(eqs_original,eqs_derivation,x_init,err_tol,Niteration)
tic;
x=x_init;
deri_inv = inv(eqs_derivation(x));
for iiIteration=1:Niteration
    x_last = x;
    %x = x-deri_inv*eqs_original(x);
    x = x-(eqs_original(x)'/eqs_derivation(x_init))';
    err_var = norm(x-x_last);
    if err_var < err_tol
        break;
    end
end
time_cost=toc;
end